Quantum interference between multiple impurities in anisotropic superconductors 
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We perform a numerical study of the quantum interference between impurities in d-wave super- 
conductors within a potential scattering formalism that easily applies to multiple impurities. The 
evolution of the low-energy local density of states for both magnetic and nonmagnetic point scat- 
terers is studied as a function of the spatial configuration of the impurities. Further we discuss the 
influence of a subdominant bulk superconducting order parameters on the interference pattern from 
multiple impurities. 
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The past few years have proved the importance of ex- 
perimental techniques which can directly test the wealth 
of information associated with modifications of the local 
density of states (LDOS) around impurities, grain bound- 
aries and vortices in strongly correlated electron systems. 
In particular scanning tunneling microscopy (STM) mea- 
surements have provided detailed LDOS images around 
single nonmagnetic 1,2 (Zn) and magnetic 3 (Ni) impuri- 
ties on the surface of the high temperature superconduc- 
tor Bi 2 Sr 2 CaCu0 4+ 5 (BSCCO). 

For conventional superconductors Yu and Shiba showed 
that as a result of the interaction between a magnetic 
impurity and the spin density of the conduction elec- 
trons, a bound state located around the magnetic impu- 
rity is formed inside the gap in the strong-scattering (uni- 
tary) limit. For anisotropic superconductors a number 
of authors generalized the Yu- Shiba approach to study 
the LDOS around single impurities 5 . It was found, for 
instance, that for a single nonmagnetic impurity there 
is only a virtual bound (or resonant) state due to the 
existence of the low-energy nodal quasi-particles. The 
one-impurity problem was recently reviewed by several 
authors 6,7 . 

Recently Hoffman et ai. 8 measured the energy depen- 
dence of the Fourier transformed LDOS images on the 
surface of optimally doped BSCCO below T c . The dis- 
persive features were explained from the point of elastic 
quasi-particle interference resulting from a single weak, 
nonmagnetic impurity 9 . This gives credence that a scat- 
tering potential picture can yield valuable predictions in 
the superconducting state of these materials. Evidence 
for quantum interference between strong scatterers has 
been observed in the CuO chains of YBa2Cu306+ x by 
Derro et ai. 10 Future experimental ability to control the 
position of the impurities on the surface of a supercon- 
ductor and perform detailed STM measurements around 
multiple impurity configurations motivates further theo- 
retical studies of the resulting quantum interference. 
Previous calculations have studied the formation of bond- 
ing and antibonding states around two magnetic impu- 
rities in s-wave superconductors 11,12 . For d-wave su- 
perconductors it was found that the interference ef- 
fects between two nonmagnetic unitary scatterers de- 



pends crucially on the distance and orientation of the 
impurities 6,7,13 . The orientational dependence arises 
from the anisotropic gap function, and provides an al- 
ternative method to identify the symmetry of the su- 
perconducting gap. Several authors 5,14 have previously 
suggested similar ideas in the case of one impurity. 
In this paper we study multiple impurity effects by ex- 
actly inverting the Gorkov-Dyson equation. In partic- 
ular, we discuss the effect of quantum interference be- 
tween: 1) nonmagnetic impurities in the strong scatter- 
ing limit, 2) nonmagnetic impurities in the case of in- 
duced subdominant superconducting order parameters, 
and 3) magnetic and nonmagnetic impurities. All the 
calculations are performed within a quasi-particle scat- 
tering framework with classical impurities 4,14 . The main 
purpose is to use quantum interference to obtain results 
that pose further tests on this approach and to illustrate 
the strong sensitivity of the LDOS on the positions of the 
impurities. 

The Greens function uj) for the unperturbed d- 

wave superconductor is given in Nambu space by 



6W(k,w) = [iw n fb-^(k)fs-A(k)f+-Atk)f_] 



(1) 



where f u denotes the Pauli matrices in Nambu space, 
To being the 2x2 identity matrix and r± = f j, ± 
T2- For a system with d x 2_ y 2-wave pairing A(k) = 
(cos(k x ) — cos(k y )). Below, A = 25meV and the 
lattice constant ao is set to unity. In this article we 
use a normal state quasi-particle energy £(k) relevant for 
BSCCO around optimal doping (14%) 

£(k) = —2t (cos(k x )+cos(k y )) — 4:t' cos(k x ) cos(k y ) — fj, (2) 

with t = 300meV, t' = -OAt and fi = -1.18*. Here t 
(t f ) refers to the nearest (next-nearest) neighbor hopping 
integral and /i is the chemical potential. 
We model the presence of scalar and magnetic impurities 
in the system by the following <5-function interactions 



[( 



(Vf+V^cl^ + iVf-V™)^ . (3) 



Here i denotes the set of lattice sites containing magnetic 
and/or scalar impurities and V^ M {Vf) is the strength of 
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the corresponding effective potential. We consider only 
the ^-component of the magnetic impurity interaction 
and ignore spin-flip scattering. 

For a single nonmagnetic impurity at the origin it is well- 
known that the scattering modifies the Greens function 

by 



6Gu(r, iuj n )=- 



G^(r,iu n ) 



G$ (r, iu v 



(4) 

Here r is the distance to the origin and G a p refers to the 
a/Jth entry of the 2x2 Nambu subspace. 
Naturally one can derive equivalent expressions for the 
LDOS modulations around several impurities. However, 
for a numerical study of the evolution of the LDOS 
for multiple impurities positioned at arbitrary lattice 
sites, we find it is easier to invert directly the real- 
space Gorkov-Dyson equation. The full Greens function 
G(r, ui) is then obtained by solving the equation 



(5) 



where the double lines indicate that the elements of this 
equation are matrices written in real- and Nambu space. 
The size of these matrices depends on the number of im- 
purities and the dimension of the Nambu space. We have 
previously utilized this method to study the electronic 
structure around vortices that operate as pinning cen- 
ters of surrounding stripes 15 . We perform the 2D Fourier 
transform of the clean Greens function G^ (k, u>) numeri- 
cally by dividing the first Brillioun zone into a 1400 x 1400 
lattice and introducing a quasi-particle energy broaden- 
ing of S = 0.5meV with S defined from iu> n — > u> + id. 
The differential tunneling conductance is proportional to 
the LDOS p(r, ui) which in turn is determined from 

p(r,u) = --Im [Gii(r,w) + G 22 (r, -u)] . (6) 

7T 

In the following we model the nonmagnetic unitary scat- 
terers with a potential V s — 700meV which gives rise to 
resonance energies around ±1.5meV in agreement with 
experiment 2 . (This is seen from the holelike resonance 
evident in the bottom LDOS scan in Fig. 5a. For a sin- 
gle impurity only one of the two resonances evident from 
Eqn. 4 has weight on the impurity site since the anoma- 
lous part of the Greens function, G\2(r,u)), vanishes at 
r = due to the symmetry of the d-wave gap.) 
For interference between two nonmagnetic unitary impu- 
rities Morr et al. 6 found strong variations in the LDOS as 
the distance between the impurities R is varied along one 
of the crystal axes. The single- impurity spectrum was ob- 
tained when R exceeds approximately 10ao- However, as 
expected for a d a .2_ J/ 2-wave superconductor, this length 
scale is much larger along the nodal directions. This is 
seen from Fig. 1. Here the density of states is measured 



above one of the impurities fixed at the origin while the 
other is moved away along the nodal (a) or anti-nodal (b) 
direction. The single impurity LDOS is obtained for R 
well above 100a . Thus only for impurity concentrations 
below 0.1% does the LDOS correspond to the expected 
result from a single strong nonmagnetic impurity. For 
weaker scatterers the decay length will be considerably 
reduced. 
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FIG. 1. DOS at (0,0) and at the single-impurity resonance 
energy ±1.5meV as a function of distance between the two 
nonmagnetic impurities separated along the (a) nodal direc- 
tion and (b) anti-nodal direction. The y-axis scale is identical 
for the two figures. 

For quantum interference between multiple fixed non- 
magnetic unitary scatterers Fig. 2 shows the LDOS as 
the STM tip is scanned along a crystal axis on which the 
impurities are positioned. 
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FIG. 2. Low bias STM scans along the horizontal axis in 
steps of ao/5 from (0, 0) (top line) to (2, 0) (bottom line). The 
scans are off-set for clarity. In (a) and (c) there are two non- 
magnetic impurities fixed at: (a) (0,0) and (1,0); (c) (0,0) 
and (2,0). In (b) and (d) there are three nonmagnetic impu- 
ritiesat: (b) (-1,0), (0,0), (1,0); (d) (-2,0), (0,0), (2,0). 

In general the resonances are split by the proximity of 
other impurities and the number of resonances is di- 
rectly proportional to the number of interfering impu- 
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rities. However, locally the density of states may be 
strongly influenced by destructive interference. For in- 
stance, for the two impurities (Fig. 2a, c) sharp reso- 
nances exist only when R = 2a as is evident from Fig. 
2c. When a third impurity is added at (—2,0) these res- 
onances appear to broaden and shift to higher energies, 
Fig. 2d. Contrary to this, Fig. 2a, b show that the ad- 
dition of a third impurity has only a minor effect when 
R = ao- 

The case of three nonmagnetic impurities is studied fur- 
ther in Fig. 3 which shows the evolution of the LDOS 
at (0, 0) as a function of the distance to a third impu- 
rity along the nodal (b,d) and anti-nodal (a,c) direction. 
The case without the third impurity corresponds to the 
topmost LDOS in Fig. 2a and 2c. 



one may utilize the quantum interference between sev- 
eral impurities as an alternative method to detect a small 
subdominant order parameter. 

For instance, tuning through a quantum phase transition 
from a d x 2_ y 2 to a d^^-Hd^ superconductor at a crit- 
ical doping level 16 , magnetic impurity concentration 17 or 
magnetic field strength 18 , a small d xy order could qual- 
itatively alter the interference pattern. For A xy (li) = 
A° xy sm(k x )sm(k y ) with A° y =5.0meV, we compare in 
Fig. 4 the LDOS at physically realizable impurity po- 
sitions to the case with A 6 =0. 
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FIG. 3. LDOS at (0, 0) as a function of distance to a 
third impurity along the antinodal (a,c) or nodal (b,d) di- 
rection. The two fixed impurities are positioned at: In (a) 
and (b): (-1,0) and (0,0); in (c) and (d): (-2,0) and (0,0). 
In (a),(c) the third impurity is positioned at (top to bot- 
tom) (1, 0),(2, 0),...,(10, 0).ln (b),(d) the third impurity is po- 
sitioned at (top to bottom) (1, 1),(2, 2),.. .,(10, 10). 

As in the case of two nonmagnetic impurities 6 ' 7 there are 
very strong variations in the final LDOS; the number of 
apparent resonances, their width and resonance energy 
depends crucially on the positions of the three impuri- 
ties. The small modulations added by the third impurity 
seen in Fig. 3a, b agree with the destructive interference 
evident from the corresponding cases seen in Fig. 2a, b. 
Contrary to this, large modulations are again seen when 
increasing the distance between the two fixed impurities 
by a single lattice constant, Fig. 3c, d. 
Recently Zhu et al. 7 suggested a careful study of the 
two-impurity problem to extract information of the bulk 
Greens function of the clean system. In the following 
we show how the quantum interference between unitary 
scatterers is strongly affected by the induction of a small 
subdominant superconducting order parameter. Thus 



Energy (me\ 



FIG. 4. Top row: DOS at (0, 0) for two nonmagnetic im- 
purities at (0,0) and (2,4). Bottom row: DOS at (0,0) (solid 
line) and (1, 1) (dashed line) for three nonmagnetic impurities 
at (-1,1), (1,-1) and (-1,-1). Pairing symmetry: (a) and (d) 
d x 2_ y 2, (b) and (e) d a .2_ !/ 2+id I y, (c) and (f) d a .2_ !/ 2+is. 

Also we show the difference between d + id and d + is 
pairing symmetry for these impurity configurations. For 
most spatial configurations the secondary pairing (id or 
is) leads to a sharpening of the resonances but at partic- 
ular positions there is a qualitative difference as shown in 
Fig. 4. For instance, the induction of d+id pairing (Fig. 
4b) can result in three apparent resonances contrary to 
the ground state with pure d^.^-wave pairing (Fig. 
4a). Similarly, by comparing the LDOS at (1,1) (dashed 
lines) in Fig. 4d-f, it is evident that the interfering scat- 
terers can provide a clear distinction between d+id and 
d+is pairing. Information of the induction of local order 
around the impurities can also be inferred from STM 
measurements of specific impurity configurations 17 . 
We turn now briefly to the study of the classical mag- 
netic impurities in d^^-wave superconductors. The 
interference of two magnetic scatterers in a s-wave su- 
perconductor was studied recently by Flatte et al. 11 For 
comparison to the nonmagnetic case we use a magnetic 
potential strength \V M \ — 700meV which does, how- 
ever, not model all magnetic impurities (e.g. Ni) on the 
surface of BSCCO 3 ' 21 . Future experiments will reveal 
whether the scattering potential formalism utilized here 
is appropriate or whether more correlated effects are 
required 19,20 . 

Fig. 5 shows the quantum interference between two uni- 



tary scatterers: (a,d) one magnetic and one nonmagnetic, 
(b,e) two parallel magnetic, and (c,f) two antiparallal 
magnetic. In all figures one impurity is fixed below the 
STM tip at the origin (0, 0) while the other is displaced 
along the: (a-c) horizontal crystal axis or (d-f) along 
the nodal direction. In Fig. 5a, d it is the nonmagnetic 
impurity that is fixed at the origin. Again the number of 
resonances, their position, amplitude and width depends 
in detail on the distance and nature of the two scatterers. 
Further, the spatial evolution of the LDOS is quite simi- 
lar for case (a),(b), and (d),(e). These are, however, very 
different from the interference between two antiparallel 
magnetic impurities (c,f) which is dominated by strong 
destructive interference at small separations along the 
antinodal direction and a surprisingly fast recovery to 
the single impurity case along both the nodal and antin- 
odal directions. 



on the bonds around the impurity site does not produce 
any qualitative changes 9 . However, on a phenomcno- 
logical level the gap suppression could allow for a com- 
peting magnetic order parameter to develop around the 
impurity. Thus, the gap suppression might be impor- 
tant for explaining the formation of magnetic moments 
around nonmagnetic impurities as seen by NMR exper- 
iments. These issues are currently controversial but the 
vast amount of information inferred from the quantum 
interference between multiple impurities may help settle 
this, and more importantly settle the validity of the scat- 
tering potential scenario versus more correlated models. 
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FIG. 5. DOS at (0, 0) for: (a,d) one magnetic and one non- 
magnetic impurity, (b,e) two magnetic (Vi' 1 = F 2 M ), (c,f) two 
magnetic (Vi M = — Vi' 1 )- The topmost graph in each figure is 
the DOS when the two scatterers are both positioned at the 
origin whereas in the lowermost at (0,0) and (10,0). Antin- 
odal separation: (a-c), and nodal separation: (d-f). 

The results presented above remain qualitative since fits 
to a specific experiment in addition to details from the 
tunneling matrix elements could also include modified 
hopping integrals around the impurities, gap suppression 
and possibly both magnetic and nonmagnetic scattering 
potentials 21 . We have checked that a gap suppression 
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